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SECTION  I 


INTRODUCTION 


The  bounded-wave  electromagnetic-pulse  (EMP)  simulator  makes  use  of 
two  parallel  finite-width  plates  as  the  guiding  structure  for  the  simulated 
EMP  (figure  1) .  One  reason  for  employing  two  parallel  plates  is  that  they 
support  a  transverse-electromagnetic  (TEM)  mode.  Another  reason  is  that 
over  a  significant  portion  of  the  region  between  the  plates,  the  TEM-mode 
fields  provide  a  good  approximation  to  the  free-space,  plane-wave  fields. 
Unfortunately,  such  a  structure  can  also  support  higher-order  transverse- 
magnetic  (TM)  and  transverse-electric  (TE)  modes  and  a  continuous  spectrum 
(refs.  1,2,3,  and  4).  The  TEM  mode  alone  is  not  sufficient  to  completely 
describe  the  total  simulator  field. 

The  properties  of  the  TEM  mode  on  two  parallel  plates  have  been 
investigated  extensively  by  the  method  of  conformal  mapping  (refs.  5,6,  and  7), 
whereas  the  higher-order  modes  and  the  continuous  spectrum  have  been  inves¬ 
tigated  only  in  some  limiting  cases.  In  reference  2,  integral  equations 
for  the  higher-order  modes  are  formulated  by  using  Green's  theorem.  The 
integral  equations  are  analytically  solved  under  the  condition  that  the 
separation  of  the  plates  is  much  larger  than  their  width  (i.e.,  narrow 
plates).  In  reference  3,  alternative  integral  equations  for  the  higher- 
order  modes  are  formulated  by  employing  Laplace  transforms  and  the  Wiener- 
Hopf  technique,  and  are  solved  for  the  plates  with  small  separation-to-width 
ratios  (i.e.,  wide  plates).  The  integral  equations  derived  in  reference  2 
are  most  useful  for  numerical  treatment  when  the  separation  of  the  plates  is 
comparable  to  or  larger  than  their  width,  whereas  those  derived  in  reference  3 
are  most  useful  when  the  separation  of  the  plates  is  comparable  to  or  smaller 
than  their  width.  The  plate  geometries  discussed  in  this  report  have 
separation-to-wldth  ratios  of  one,  two  and  three.  The  integral  equations 
derived  in  reference  2  are  thus  more  appropriate.  In  this  report,  numerical 
results  for  the  propagation  constants  and  field  distributions  will  be  given 
for  the  TM  modes.  The  TE  modes,  which  are  more  highly  attenuated  away  from 
the  launching  region  (ref.  2),  will  be  discussed  in  an  appendix. 


There  is  not  much  information  available  regarding  the  continuous 
spectrum  of  the  two-parallel-plate  simulator.  In  reference  1,  an 
asymptotic  analysis  has  been  given  to  calculate  the  contribution  of  the 
continuous  spectrum  to  the  total  field  of  two  Infinitely  long  parallel 
wires,  which  may  be  considered  as  a  limiting  case  of  a  two-parallel-plate 
simulator.  In  this  report,  a  preliminary  asymptotic  estimation  of  the 
continuous  spectrum  contribution  to  the  TM  field  at  a  fixed  frequency  will 
be  given. 

To  solve  the  integral  equations  derived  in  reference  2,  one  first 
transforms  the  integral  equations  into  the  Fredholm  integral  equations  of 
the  second  kind  by  using  Carleman's  formula  for  singular  integral  equations 
(ref.  8).  The  resulting  integral  equations  are  transformed  further  into 
matrix  equations  which  can  be  solved  numerically  by  expanding  the  unknown 
functions  in  terms  of  Chebyshev  polynomials.  From  the  numerical  solutions 
of  the  matrix  equations,  the  propagation  constants  and  field  distributions 
of  the  higher-order  modes  as  well  as  the  properties  of  the  continuous 
spectrum  can  be  obtained  by  some  straightforward  calculations. 

It  should  be  mentioned  that  although  the  results  in  this  report  are 
obtained  for  the  infinitely  long  plates,  they  can  be  directly  applied  to  the 
real  simulators  where  the  lengths  are  finite. 


SECTION  II 


INTEGRAL  EQUATION  FORMULATION 

Two  infinitely  long,  perfectly  conducting,  parallel  plates  of  finite 
width  are  shown  in  figure  2.  The  width  of  each  plate  is  2w  and  the  distance 
separating  the  plates  is  2h.  A  coordinate  system  is  introduced  such  that 
the  z-axis  coincides  with  the  axis  of  the  structure  and  the  x-y  plane  is 
the  transverse  plane,  with  the  x-axis  parallel  to  the  plates. 

The  transverse  field  components  E^x,y,£,s)  and  H^x.y.c.s)  in  the 

Laplace  transform  domain  (i.e.,  s,  £  domain)  are  related  to  the  longitudinal 

field  components,  E  (x,y,£,s)  and  H  (x,y,£,s),  via  (ref.  2) 
z  z 

—2  -2  * 

E  (x,y,£,s)  -  -  £p“  7  E  (x,y,£,s)  -  sp  p"  zx  7  H  (x,y,£,s) 

C  Z  Z  O  t  Z 

^(x.y.C.s)  *  -  Cp“2VtHz(x,y,;,s)  +  seQp“2  z  x  V^U.y.C.s) 

2  2  2  2 

where  the  factor  exp(£z  +  st)  has  been  suppressed,  p  » s  /c  -  £.  ,  c  is  the 
vacuum  speed  of  light  and  Vt  is  the  gradient  in  the  transverse  direction. 

From  equation  1,  it  is  obvious  that  the  fields  can  always  be  decomposed 

into  two  parts,  the  TM  fields  with  H  *0  and  the  TE  fields  with  E  *0.  Each 

z  z 

part  will  be  discussed  separately  in  the  sequel. 

1.  TM  FIELDS 

As  has  been  discussed  above,  a  knowledge  of  Ez  is  sufficient  for  the 

determination  of  the  TM  field  distributions.  From  reference  2,  E  (x,y,C,s) 

z 

is  given  by 

Ez(x,y,£,s)  »  |  G(x,y,x' ,h;p)f (x' ,h)dx'  +  j”  G(x,y,x' ,-h;p)f (xr ,-h)dx' 

-w  -w 

G(x,y,x'  ,y'  ;p)  -  ±  kJP/(x-x’)2  +  (y-y')2  )  (2)  j 

f  (x,±h)  ■  lim  (■!—  E  (x,±h+e,£,s)  -  t—  E  (x,±h -  e,£,s))  1 


where  Kq  is  the  modified  Bessel  function  of  the  second  kind  and  f(x,±h) 
satisfy  the  following  set  of  integral  equations 

fW  fW 

G(x,h,x' ,h;p)f (x' ,h)dx'  +  G(x,h,x' ,-h;p)f (x' ,-h)dx'  *  a(x,h) 
-w  -w 

rrw 

G(x,-h,x' ,h;p)f (x' ,h)dx'  +  j  G(x,-h,x' ,-h;p)f (x',-h)dx'  »a(x,-h) 


for  |x|  <.  w,  with  a(x,±h)  »  -  (x,±h,£,s)  being  the  incident  longitudinal 

electric  fields  at  the  plates. 

2.  TE  FIELDS 

The  TE  field  distributions  can  be  derived  completely  from  Hz(x,y,c;,s) 
which  is  given  by  the  following  formula  (ref.  2) 

spoH2(x,y,C,s)  *■  -  j  ||p-  (x.y.x'.y'jpJgCx*  ,h)dx' 

-w  *-  y  -V '  =  h 


-  j  [§p-  (x,y ,x',y '  ;p)g(x' ,-h)dx'j 


y'  -  -h 


where  g(x,±h)  are  defined  by 


(suQ)  1g(x,±h)  -  lim  (  H  (x,±h  +  e,C,s)  -  H  (x,+h  -  e,£,s)) 


and  satisfy  the  following  set  of  differential-integral  equations 

G(x,h,x,,h;p)g(x',h)dx'  +  G(x,h,x' ,-h;p)g(x' ,-h)dx’^ 


B(x,h)  ,  for  |x|  <_  w 


j-p2^!  G(x,-h,x' ,h;p)g(x' ,h)dx'  +  G(x,-h,x'  ,-h;p)g(x'  ,-h)dx^ 


with  8(x,±h)  *  p^E^nc(x,±h,5,s)  +  5 3Einc (x,±h,? ,s) /3x  being  the  source  terms. 

X  2 

Although  equations  3  and  5  look  extremely  complicated,  they  can  be 
simplified  by  observing  that  for  most  parallel-plate  simulators  the  source 
terms  on  the  right-hand  sides  of  equations  3  and  5  satisfy  the  following 
conditions 


ct(x,h)  ■  -  a(x,-h)  »  a(-x,h)  =  a  (x) 
8(x,h)  =  -  g(x,-h)  =*  -8(-x,h)  =  S  °(x) 


(6) 


Accordingly ,  one  has 

f(x,h)  =  -  f (x,-h)  =  f(-x,h)  =  f-e(x) 
for  the  TM  fields,  and 

g(x,h)  =>  -  g(x,-h)  *  -  g(-x,h)  *  g"*°(x) 

—3  — o  —0  —  o 

for  the  TE  fields.  The  superscript  "  -  "  in  a  ,8  ,  f  ,  g  is  used  to 

indicate  that  all  these  functions  a,  8,  f,  g  are  antisymmetric  with  respect 
to  y,  while  the  superscript  "e"  or  "o"  is  used  to  indicate  that  the  functions 
are  either  even  or  odd  functions  of  x.  Under  the  above  conditions  of  equation  6, 
the  two  equations  of  either  integral  equation  set  3  or  5  become  identical. 

In  the  following  sections,  the  simplified  equations  will  be  used  to 
investigate  the  properties  of  the  higher-order  modes  and  the  continuous 
spectrum  of  the  two-parallel-plate  guiding  structure. 
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SECTION  III 


MATRIX-EQUATION  FORMULATION  FOR  THE  TM  FIELDS 

la  this  section,  the  integral  equations  given  by  equation  3  for  the  TM 
fields  will  be  transformed  into  matrix  equations  by  expanding  the  unknown 
functions  in  terms  of  Chebyshev  polynomials. 

On  account  of  the  source  condition  of  equation  6  the  coupled  integral 
equations  3  are  reduced  to  a  single  integral  equation 


fo(x,h,x' ,h;p)  -  G(x,h,x' ,-h;p)  j  f  e(x')dx'  *  a  ®(x),  |x|  £w  (7) 

-w  ' 

It  is  easy  to  see  that  the  kernel  of  the  integral  equation  7  has  a  loga¬ 
rithmic  singularity.  After  separating  out  the  singular  term  and  normalizing 
the  variables  in  the  following  manner 


x  =  w£,  x'  «  w 
p  =  y/w,  h  *  w  H 


(8) 


the  integral  equation  7  becomes 

£n( | C-C  |)f-e(£')d£' 

-1  x 

=*  2 it  f  (m(5-5';y)  -N(5-C,;Y))f-e(S,)d£’  -  27Ta"e(C),  for  |c|  <  1 
*  1 


where 

m(c-£';y)  -  ^  ko(y|c-c'|)  Ende-e’l) 

N(C-C’;y)  -  Ko  (y/c-5')2  +  4h2  ) 

Here,  both  M  and  N  are  regular  functions  of  £-  C'* 


10 


■  -..*3 


To  facilitate  the  numerical  solution  of  the  Integral  equation  9,  it 
is  desirable  to  transform  it  into  a  Fredholm  integral  equation  of  the 
second  kind.  To  do  this,  Carleman's  formula  for  integral  equations  with 
logarithmic  kernels  is  used  to  obtain  (ref.  8) 

f"e<5>  +  [^KCC.C'jy)  -  LU.C’jy))  f"ea*)dC'  -  S_e(g),  |C|  <1  (10) 


where 


k(C,C’;y) 


— J- f1  mizJLnl  dV, 

it  A  -  g2  In  2  -1 

M1  (g"  -  g  *  ;y)  -  g'*^ 


\\ 


tt/I  -  g*  -1 


l  -  C 


dg" 


L(C,g' ;y)  - 


r  A  -  g2  in  2  _1  Jl  ~  F" 


f  dg” 


tt/I  -  g 


j1  N,q,,-g,;Y)/i-g'2  ,c.. 
U  5-g" 


2  J-1 


Here,  £  denotes  the  principal-value  integral  and  the  prime  in  M’ ,  N' ,  a"6' 
denotes  differentiation  with  respect  to  the  first  argument. 

The  integral  equation  10  can  now  be  reduced  to  a  set  of  algebraic 
equations  suitable  for  numerical  computation.  The  mathematical  properties 


of  K(5,5';y),  L(£,5';y)  and  S-e(C)  together  with  the  edge  conditions  and 
symmetry  properties  of  f  (?)  suggest  the  expansion 


f‘*(5) 


'1-5 


2  n=0 


l  Cvo 


(id 


where  T.  (?)  are  Chebyshev  polynomials  of  the  first  kind.  Then,  with  the  aid 
zn 

of  the  orthogonality  of  the  Chebyshev  polynomials,  the  integral  equation  10 
is  transformed  to  the  following  set  of  algebraic  equations 


where 


_e  +  l  ( Ke  -  L®  ) 
n  Ln  \  nm  nm  / 
m*u 


Ke  -  L®  )  f"e  -  S_e 
nm  nm  /  m  n 


nm  ire 


f. 


I  - — - 2“ - dC’dC 


n 


A- e'2 


(12) 


.  o  r1  ! 1  ;y)t?  (5)T~  (5') 

Le  a  —2—  j  - 25 - 2m -  dc,d? 

“  ff6n  U  U 


,-e  2 

5  *  - 

n  ire 


fi 


s"e(5)T2n(5)d5 


and 


'2,  n *  0 
'1,  n^O 


Or,  in  matrix  form 


[«  +  K®  -  L® 

L  nm  nm  nm 


][cHc] 


(13) 


where  6  is  the  Kronecker  delta,  which  is  zero  when  n  is  not  equal  to  m 
nm 

and  unity  when  n  equals  m,  and  [r  ]  Is  a  matrix  whose  elements  are  rnm. 


0  g 

Both  K  and  L  are  complicated  triple  integrals.  They  can  be 

nm  nm 

simplified  by  first  integrating  over  £  and  then  by  making  the  change  of 
variables  -  2q  and  £"+£'  *  2n’,  to  get 


TTST  f  M(2n;Y)G*(n)dn  n-0 

o 

f  M' (2n;y)F®  (n)dn  n  >  1 

tr  j  nm 

'o 

TtoT  |  N(2n;Y)G*(n)dr1  n-0 

0 

(  N'(2n;y)F®  (n)dn  n  >  1 

it  i  nm 

J0 

where 


(14) 


(15) 


and  Un(n)  is  the  Chebyshev  polynomial  of  the  second  kind.  By  using  the  expli¬ 
cit  formulas  for  the  Chebyshev  polynomials  in  the  integrals  of  equation  15, 

6  6 

G  (n)  and  F  (n)  are  further  transformed  into  sums  of  complete  elliptic 

tn  nm 

integrals  which  are  more  suitable  for  numerical  computation,  viz., 


G*(n)  -  l  l  t*  b.(0,2k)n2k'2jD,(n) 
k-0  j-0  J  3 

ra  n  k+A-1  .  7 

F®  (n)  ■  u*  t®  l  (b . (21-1,  2k)  -  rTb.(2Z+l,  2k)  x  (16) 

0111  k-0  Jl-1  a*  J-0  V  j  j  ’ 

*  n2i  +  2k-  2j  -  1^  (TJ)  _  n  bfc+Jl(2i+l,  2k)Dk+£(n) 


13 


where 


e  .  ..n-A  (n  +  A  -  1)  !  „2A  -  1 

\t  '  <-l)  <0 - 1) i72t - 'i) i  2 


te  -  (-d®-*  te.-1-*-.1?*®  22k  te  - 1 

Cmk  '  A;  (b  -  k) !  (2k) !  *  Coo  A 


b^k.A) 


r  ,  nA+i-2j  _ klJJ _ 

i,  11  t!  (k  -  i) !  (2j  -  i) !  (A  +  i  -  2j ) ! 


-  max{0,2j  -  a}  i2  -  min{k,2j} 

and  Che  funeciona  Dj(n)  are  determined  from  che  following  recursion  formulae 

2 

Dj+2(n)  "  4  2jT3  ^ 1  +  r>2)Dj+1(n)  -  (l-n2) 

D1(n)  -  2(1  +  n) |f( (1  -  n)/(l  +  n))  -  E((l-n>/(l  +  n))j  (17) 

Do(n)  -  2F  ((l  -  n)/(l  +  n))/ (1  +  n) 

Here,  F(n)  and  E(n)  are  complete  elliptic  Integrals 

F(n)  »  f  (l- n2sin2^)'1*d<> 

E(rj)  -  |  ( 1  -  n2sin2<^)Jsd^ 

0 

The  solution  of  the  matrix  equation  13  is  simply  given  by 

[f"*l  -  \s  +  K*  -  L®  Hs-1  (18) 

L  n  J  L  o®  nm  nmj  L  n  J 

Insertion  of  this  solution  into  equations  11,  2  and  1  gives  the  Laplace 
transform  domain  TM  fields. 


14 


In  the  next  two  ■actions  certain  important  properties  of  the  TM  fields 
will  be  discussed  by  studying  the  singularities  of  equation  2  in  the 
complex  (-plane  at  a  fixed  ■  ■  jut. 


SECTION  IV 


PROPAGATION  CONSTANTS  AND  FIELD  DISTRIBUTIONS  OF  DISCRETE  TM  MODES 


In  Che  previous  sections,  Che  formulas  required  Co  calculate  Che  TM 
fields  of  two  parallel-places  are  obtained.  At  a  given  frequency  ui,  the 
field  distributions  are  calculated  from  the  Inverse  Laplace  transform 
Integrals 


1  E(x,y,z,jw) 

E(x,y ,c»jw) 

j  H(x,y,z, jw) 

H(x,y,;,ju) 

(19) 


where  C^  Is  the  path  of  Integration  In  the  complex  £ -plane  shown  In  figure  3 
and  E(x,y,C.jw)  and  H(x,y ,C,jw)  are  given  by  equations  1  and  2.  An  examina¬ 
tion  of  equation  2  shows  that  In  the  complex  C-plane,  there  are  branch  points 

at  ;  ■  i  liu  and  poles  at  C,  „  where  det  6  +  Ke  -  Le  “0.  For  the 

k.  >  V*  run  tun  run  j 

branch  cuts  shown  In  figure  3,  the  contour  C^  can  be  deformed  to  the  left 
half-plane  for  field  points  at  z  >  0  and  sources  at  z  <.0.  The  contour 
Integral  along  C  Is  thus  reduced  to  the  Integral  along  the  branch  cut  B_ 

(the  so-called  continuous  spectrum  contribution)  plus  the  residues  at  the 
poles  ^ (the  so-called  modal  fields).  The  reason  why  two  Indices  are 
assigned  to  ^  ^  will  become  clear  later  when  numerical  results  are  obtained. 

Each  modal  field  Is  required  to  be  outgoing  in  the  transverse  direction 
and  decaying  In  the  +z  direction.  Thus,  in  the  branch  shown  in  figure  3,  the 
poles  can  exist  only  in  the  region  where  -u>/c  <  Im(£)  <  0  and  Re(c)  <  0. 

It  is  also  observed  that  p^  ^  have  negative  real  parts,  so  that  the  modal 
field  distributions  increase  Indefinitely  in  the  transverse  direction. 

Hence,  the  branch  shown  in  figure  3  is  appropriate  only  when  the  field 
points  are  close  to  the  plates.  In  the  remaining  part  of  this  section,  the 
propagation  constants  and  field  distributions  of  the  TM  modes  will  be 
discussed,  while  an  estimation  of  the  continuous-spectrum  contribution  will 
be  relegated  to  the  next  section. 

To  calculate  the  propagation  constants  of  the  TM  modes,  one  has  to 

find  first  p,  .  from  the  equation  det  5  +  Ke  -  Le  *0.  With  these 

k, i  n  L  nm  run  nm  J 
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Figure  3.  The  Branch  of  p 


-  j/c2  +  <*> 


2/2 


/c  .  In  this  Branch,  Re(p) 


for  {5:  Re(c)  <  0  and  -u»/c  <  Im(c)  <  0}  and  {c:  Re(?)  > 


Pk  ^  values  (the  so-called  transverse  propagation  constants),  the  corresponding 
longitudinal  propagation  constants  ?k  ^  are  simply  given  by 

1  /<4,*  +  “2/cJ  <20> 

When  H  •  h/w  »  1,  i.e.,  the  plates  are  extremely  narrow,  there  exist 
possible  pkys  for  which  Iy^J  -  |pk>iw|  «  1  and  det  [  5^  +  -  LjJ-0. 

Actually,  under  the  condition  |yk  ^ |  <<  1,  one  can  show  that  the  yk  ^'s 
satisfy 

det  f <5  +  K®  -  L®  1  -  det  +5  6  (K®  -  L®  )  1 

L  tun  nm  nmj  L  o®  no  mo  oo  ooj 

*  1  +  (in  2-  in(rYkjl>  -  V2\,JlH)  )  /ln2  “  0  (21) 


where  r  ■  1.781  is  the  exponent  of  Euler's  constant.  Equation  21  is  the 
same  as  that  obtained  in  reference  2,  where  its  solutions  are  also  given. 

When  the  width  of  the  plates  is  comparable  to  their  separation,  which 

is  the  case  of  interest  here  where  H  -  1,2,3,  one  must  resort  to  numerical 

methods  to  solve  the  equation  det  <5  +  K  -  L®  •  0  for  the  transverse 

n  L  nm  nm  nmj 

propagation  constants  p,  . .  The  method  selected  is  first  to  locate  the 

ICf 

approximate  positions  of  the  zeros  of  the  determinant  from  the  constant- 
magnitude  contours  in  the  complex  p-plane.  Then  these  approximate  positions 
are  used  as  the  starting  points  for  the  Newton-Raphson  method,  used  to 
search  for  more  accurate  solutions  for  the  transverse  propagation  constants 
£.  It  is  found  that  three  terms  in  the  expansion  11  for  f  (5)  are 
sufficient  to  obtain  resonably  accurate  transverse  propagation  constants 
of  the  first  several  modes  for  H  ^  1.  For  each  transverse  propagation 
constant  pk  the  fQ  are  then  determined  within  a  multiplicative  constant 
from  the  homogeneous  equation  of  equation  13.  These  f  values  can  in  turn 
be  used  in  equations  1,  2  and  11  to  calculate  the  corresponding  modal  field 
distributions. 


The  transverse  propagation  constants  p,  .  of  the  first  several  TM 
modes  for  H  -  1,2,3  calculated  by  the  above  described  approach  with  the 
aid  of  a  CDC  6600/7600  computer  are  tabulated  in  table  1.  The  ^ 
values  as  functions  of  u>  for  the  lowest  TM  modes  with  k*0,l  and  1*1,2 
are  plotted  in  figures  4,  5  and  6.  Four  curves  of  p^  i  m  1,2, 3, 4  are 
also  presented  in  figure  7  for  H-values  ranging  from  50  to  about  2.  The 
numerical  results  of  the  modal  field  distributions  for  the  TM^  ,  modes 
with  k*0,l,  )l-l,2  are  given  in  figures  8  through  19.  The  field  distribu¬ 
tions  are  plotted  in  terms  of  constant-value  contours  of  the  real  parts, 
imaginary  parts  and  magnitudes  of  the  normalized  field  components. 

From  the  field  distribution  plots  of  the  TM^  ^  modes,  especially  the 
modes  with  l  >  k,  it  is  observed  that  in  the  working  volume  of  the  simulator 
(i.e.,  x/w,  y/h  <  1)  the  fields  vary  almost  sinusoidally  as  functions  of  x 
and  y  with  periods  of  2w/k  and  2h /l  respectively.  The  indices  k,&  used  in 
the  TM^  ^  modes  thus  characterize  the  field  variations  in  the  x,y  directions 
and  the  use  of  two  indices  is  justified.  In  this  report,  results  are  given 
only  for  k-0  and  1.  It  is  believed  that  if  one  goes  even  farther  away  from 
the  imaginary  axis,  more  TM^  ^  modes  can  be  found  for  k  >,  2.  However,  for 
the  modes  with  k  >_  2,  the  corresponding  longitudinal  propagation  constants 
?,  will  have  decay  constants  so  large  that  those  modes  become  less 

K)  M 

important . 


Table  1 

TRANSVERSE  PROPAGATION  CONSTANTS  OF  TM  MODES 


k 

-  0 

k 

*  1 

Re(p0jjlw) 

R6(p1,£W) 

Im(pMw> 

-0.1274 

3.2879 

-1.5701 

1.1930 

-0.0590 

6.3890 

-0.9827 

4.3512 

-0.0353 

9.5060 

-0.5034 

7.1986 

-0.0245 

12.6321 

-0.2935 

10.1311 

-0.1834 

1.6686 

-1.2325 

0.7914 

-0.1260 

3.2448 

-1.1963 

2.2510 

-0.0778 

4.8067 

-0.9098 

3.9028 

-0.0575 

6.3662 

-0.6898 

5.4269 

-0.1779 

1.0974 

-1.0351 

0.5888 

-0.1299 

2.1721 

-0.9925 

1.4961 

-0.1006 

3.2207 

-0.9207 

2.5827 

-0.0814 

4.2662 

-0.7896 

3.6733 

*  Atmans* 


Figure  4.  Frequency  Variation  of  Longitudinal  Propagation  Constants 
5.  of  Higher-Order  TM  Modes  for  k  =  0,l,  1  =  1,2  when  h/w 

•V}  JC 
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Figure  5.  Frequency  Variation  of  Longitudinal  Propagation  Constants 

r  of  Higher-Order  TM  Modes  for  k*0,l,  £*1,2,  when 
k,  £ 

h/w  »  2. 
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Figure  7.  The  Normalized  Transverse  Propagation  Constants  p0  {w  as 
Functions  of  h/w  for  l  =  1,2, 3, 4. 


for  Negative  Values. 


x/w 


Figure  8c.  Constant  Value  Contours  for  Normalized  Field  Component  | E^/w | 
of  the  TMq  ^  Mode  when  h/w  -  1.  The  Fields  are  Normalized 
so  that  |p2/s|/|ExJ2+  |Ey|2  -  1  at  x-y-0. 
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Figure  8e.  Constant  Value  Contours  for  Normalized  Field  Component _ 

/(Im(p2Ex/0)2+  (Im(p2Ey/c))2  -  /(Im(p2Hy/seo))2  +  (InKp2^/ 

of  the  TMq  ^  Mode  when  h/w  »  1.  The  Fields  are  Normalized 

c*  1  2  .  I  p  I  2  1  ~ 


so 


Figure  9b. 


Constant  Value  Contours  for  Normalized  Field  Component 

Im(E  /w)  of  the  TMn  ,  Mode  when  h/w  ■  1.  The  Field*  are 

z  / - 5 — - y 

Normalized  so  that  |p  /c|/|E  |  +|E  |  -  1  at  x  -  y  -  0. 

*  y 

Broken  Lines  are  for  Negative  Values. 


Figure  9c. 


^“7,  — u-d 

l  *'"1  °‘  th‘  m0,2  «*■«•  *>«  h/.  -  , 
°n“U!,d  ••  th.t  ip^ci/TTCTrrr 


Field  Component 
The  Fields  are 
■  1  at  x  ■  y  »  o. 


x/w 


Figure  9e.  Constant  Value  Contours  for  Normalized  Field  Component 


/(Im(p2E  / ?))2  + (Im(p2E  / c))2 


/(Im(p2H  /se„))2  +  (Im(p2Hv/se^))2 


of  the  TMq  ^  Mode  when  h/w  ■  ] 
that  |P2/(;Mex|2+  |Ey I2  -  1  at  x-y 


y'  o' '  '  ‘  x  o 

The  Fields  are  Normalized  so 


0. 
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Figure  10a.  Constant  Value  Contours  for  Normalized  Field  Component 

Re(E  /w)  of  the  TM.  .  Mode  when  h/w  =  1.  The  Fields  are 
2  /  2  2 

Normalized  so  that  Ip  /c|/|E  |  +  |E  |  =  1  at  x“y  =  0. 

x  y 

Broken  Lines  are  for  Negative  Values. 
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/(Im(p  Ex/C))  +(Im(p  yO)  -  /(Im(p  Hy/seo))  +  (Im(p  Hx/seo)) 
of  the  TM^  ^  Mode  when  h/w  ■  1.  The  Fields  are  Normalized  so 
that  |p2/5|/|Ex|2+  ) E  ] 2  -  1  at  x-y»0. 


I 


Figure  11a.  Constant  Value  Contours  for  Normalized  Field  Component  Re(Ez/w) 
of  the  TM,  „  Mode  when  h/w  -  1.  The  Fields  are  Normalized  so 

2  /  o  ? 

that  |p  /c|/|eJ  +  |E  I  *  1  at  x-y*0.  Broken  Lines  are 
x  y 

for  Negative  Values. 
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Figure  lib.  Constant  Value  Contours  for  Normalized  Field  Component 

Im(E  /w)  of  the  TM,  ,  Mode  when  h/w  -  1.  The  Fields  are 
z  /  2  2 

Normalized  so  that  |p2/c|/|Ex|  +  [ E  |  ■  1  at  x*y*0. 

Broken  Lines  are  for  Negative  Values. 


Figure  lid.  Constant  Value  Contours  for  Normalized  Field  Component 


/(Re(p2Ex/?))2  +  (Re(p2Ey/i;))2  =  /(Re(p2Hy/seo)  )2  +  (Re(p2Hx/seo) ) 
of  the  TIL.  „  Mode  when  h/w  ■  1.  The  Fields  are  Normalized  so 


that  I p 


+  E. 


1  at  x  -  y  *  0. 
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Constant  Value  Contours  for  Normalized  Field  Component 

Re(E  /w)  of  the  TM«  ,  Mode  when  h/w  ■  2.  The  Fields  are 

®  /  2  —  2 

Normalized  so  that  |p  /q|/|E  (  +[E  |  *  1  at  x*y*0. 

x  y 

Broken  Lines  are  for  Negative  Values . 


x/w 


Figure  12d.  Constant  Value  Contours  for  Normalized  Field  Component 

/(Re(p2E  /c))2  +  (Re(p2E  /O)2  *  /(Re(p2H  /se  ))2  +  (Re(p2H  /sc  )) 

A  y  V  U  A  U 

of  the  TMq  ^  Mode  when  h/w  =  2.  The  Fields  are  Normalized  so 

that  | p2/ S |  /[ E  | 2  +  | E  [ 2  •  1  at  x  ■  y  *  0. 
x  y 
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2 


Figure  l2e.  Constant  Value  Contours  for  Normalized  Field  Component 

/(Im(p2E  /C))2  + (Im(p2E  /S))2 =  /(Im(p2H  /se  ))2 +  (Im(p2Hx/seo) ) 
x  y  j 

of  the  TMq  .Mode  when  h/w  “2.  The  Fields  are  Normalized  so 
that  [p2/c(/|Exl2+  |Ey(2  -  1  at  x-y-0. 


Figure  12f.  Constant  Value  Contours  for  Normalized  Field  Component 


Mode  when  h/w  «*  2.  The  Fields  are  Normalized  so  that 


x/w 


Figure  13b.  Constant  Value  Contours  for  Normalized  Field  Component 

Im(E_/w)  of  the  TM-  ,  Mode  when  h/w  *  2.  The  Fields  are 
^  0  j  2-  /  2  o' 

Normalized  so  that  |p  /;|/jE  |  +  |E  |  ■  1  at  x*y*0. 

x  y 

Broken  Lines  are  for  Negative  Values. 
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int  Value  Contours  fc 

,  tMq  2  Mod^jAsaJii 

at 


(Re(p2E  /c))2  +  (Re(p2Ew/0)2 
A  y 


(Re(p2Hy/seo))2  + (ReCp^/s^))2 


of  the  TM_  ,  Mode  when  h/w  *  2.  The  Fields  are  Normalized  so 


that  |p  /cl 


1  at  x  ■  y  ■  0. 


1 


x/w 

Figure  I3e.  Constant  Value  Contours  for  Normalized  Field  Component 

/(Im(p2Ex/c))2+ (Im(p2Ey/;))2  -  /(Im(p2Hy/seo))2  +  (Im(p2Hx/seo))2 
of  the  TMq  t  Mode  when  h/w  ■  2.  The  Fields  are  Normalized  so 
that  1  p2/c  |  /j  Ex  | 2  +  |  Ey  | 2  >  1  at  x-y-0. 
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Figure  14a. 


Constant  Value  Contours  for  Normalized  Field  Component 


Re(E_/w)  of  the  IM,  ,  Mode  when  h/w  =  2.  The  Fields  are 
z  l  9  j.  2  i  2  o 

Normalized  so  that  |p  /c|/|Ex|4+|E  j  *  1  at  x»y*0. 

Broken  Lines  are  for  Negative  Values. 


Figure  14b.  Constant  Value  Contours  for  Normalized  Field  Component 


Im(Ez/w)  of  the  TM^  Mode  when  h/w  -  2.  The  Fields  are 
Normalized  so  that  |p2/c  j  /\  eJ  2  +  |Ey|2  *  1  at  x-y-0. 
Broken  Lines  are  for  Negative  Values* 


Figure  14c.  Constant  Value  Contours  for  Normalized  Field  Component  lE^/w 
of  the  ^  Mode  when  h/w  *  2.  The  Fields  are  Normalized 
so  that  |  p2/^  [  /j  E  |2+  |  E  1 2  -  1  at  x»y*0. 


x/w 

Figure  14e.  Constant  Value  Contours  for  Normalized  Field  Component 

/(Im(p2Ex/i;))2  +  (Im(p2Ey/c))2  =  /(Im(p2Hy/seo))2  +  (Im(p2Hx/seo) ) 2 
of  the  TM^  Mode  when  h/w  =  2.  The  Fields  are  Normalized  so 
that  |p2/?i/|Ex|2+  |  Ey  1 2  =  1  at  x  =  y*0. 
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0.5 


1.5 


2 


x/w 


Figure  14 f. 


Constant  Value  Contours  for  Normalized  Field  Component 


|p2/c|/|Ex|2+  |Ey|2  =  |p2(s£0)-1|/!hx|2  +  |Hy|2  of  the  TM^ 
Mode  when  h/w  °  2.  The  Fields  are  Normalized  so  that 
lp2/4|/iEx!2+  |Ey|2  3  1  at  x  =  y  =  0. 
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0  0.5  I  1.5  2 


x/w 

Figure  15c.  Constant  Value  Contours  for  Normalized  Field  Component  | Ez /w | 
of  the  TM^  2  Mode  when  h/w  *  2.  The  Fields  are  Normalized 
so  that  |p2/?|/|E  | 2  +  J E  | 2  -  1  at  x-y-0. 

x  y 


I 

r 


i' 


x/w 

Figure  15d.  Constant  Value  Contours  for  Normalized  Field  Component 

/(Re(p2Ex/i;))2+  (Re(p2Ey/5))2  -  /(Re(p2Hy/seo) )2  +  (Re(p2Hx/seo) ) 
of  the  TM^  2  Mode  when  h/w  =  2.  The  Fields  are  Normalized 
that  |p2/cl/|E  1^+  jE  |2  =■  1  at  x-y-0. 

a  y 
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so 


Figure  I5e. 


Constant  Value  Contours  for  Normalized  Field  Component 
/(lm(p2E  JO)2  +  (Im(p2E/5>)2  -  /(lm(p2H  /se ~)7  +  (Im(p 

x  y  / 


of  the  TM.  „  Mode  wi — 


TU  —  -1  J, 


learnt 
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Figure  16a.  Constant  Value  Contours  for  Normalized  Field  Component 

Re(E_/w)  of  the  TM„  ,  Mode  when  h/w  ■  3.  The  Fields  are 
2  u ,  x.  _  r~ — o  2 

Normalized  so  that  |p  /c  | /|  E  |  +  |  E  |  ■  1  at  x-y“0. 

x  y 

Broken  Lines  are  for  Negative  Values. 


x/w 

Figure  16b.  Constant  Value  Contours  for  Normalized  Field  Component 

Im(E  /w)  of  the  TMn  ,  Mode  when  h/w  *  3.  The  Fields  are 
z  U.J--  / - 2 — - 2 

Normalized  so  that  |p  /g|»|e  |  '+|E  |  *  1  at  x*y=*0. 

x  y 

Broken  Lines  are  for  Negative  Values. 


Figure  16c. 


x/w 

Constant  Value  Contours  for  Normalized  Field  Component  | Ez/w | 
of  the  TMq  ^  Mode  when  h/w  ■  3.  The  Fields  are  Normalized  so 
that  !p2/?j/|Ex|2+  |Ey|2  -  1  at  x-y-0. 


x/w 

Figure  16d.  Constant  Value  Contours  for  Normalized  Field  Component 

/(Re(p2E  /O)2  +  (Re(p2E  /O)2  -  /(Re(p2H  /st  )  )2  +  (Re(p2H /se  )  )2 
x  y  y  o  x  o 

of  the  TMq  ^  Mode  when  h/w  ■  3.  The  Fields  are  Normalized  so 
that  |p2/i!/|Ex|2+  j Ey ] ^  -  1  at  x»y-0. 
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Figure  17a.  Constant  Value  Contours  for  Normalized  Field  Component 

Re(E  /w)  of  the  TMn  _  Mode  when  h/w  =  3.  The  Fields  are 

0,2  - ~ - r 

Normalized  so  that  |p  /c|/|e  |z+|e  |  =  1  at  x  =  y  =  0. 

x  y 

Broken  Lines  are  for  Negative  Values. 


x/w 


Figure  17b.  Constant  Value  Contours  for  Normalized  Field  Component 

Im(E  /w)  of  the  TMn  ,  Mode  when  h/w  =  3.  The  Fields  are 
Normalized  so  that  j  p  /  C  | »  [  |  +|e|  =  1  at  x  =  y  =  0. 

Brrxen  Lines  are  for  Negative  Values. 
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Figure  17c.  Constant  Value  Contours  for  Normalized  Field  Component  lEz/w 
of  the  TMn  0  Mode  when  h/w  *  3.  The  Fields  are  Normalized 

U ,  ^  r— 

so  that  | p^/5  |  /  | Ex|  +  | E  |  =  1  at  x  =  y  =  0. 
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x/w 

Figure  l7f.  Constant  Value  Contours  for  Normalized  Field  Component 

!P2/S|/!EX|2+  |Ey  I2  =  jp2Cseo)“1f /|Hx|2  |Hy|2  of  the  TMQ<2 

Mode  when  h/w  *»  3.  The  Fields  are  Normalized  so  that 
lp2/?|/SEx|2+  |Ey|2  =  1  at  x«y«0. 
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Figure  18e.  Constant  Value  Contours  for  Normalized  Field  Component 

/(Im(p2E/c))2+  (Im(p2E  /c>)2-  /(Im(p2H  /se  ))2  +  (Im(p2H  /se  )) 
x  y  x  o  y  o 

of  the  TM^  ^  Mode  when  h/w  =  3.  The  Fields  are  Normalized  so 
that  |  p2/ ?  |  /|E  | 2  +  |E  | 2  -  1  at  x*y  =  0. 

x  y 
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o  I — 1  1  ■  t  L  ~  1  I 

0  0.5  I  1.5  2 

x/w 

Figure  19b.  Constant  Value  Contours  for  Normalized  Field  Component 

Im(E  /w)  of  the  TM,  ~  Mode  when  h/w  =  3.  The  Fields  are 
z  i»2  t  ^  T 

Normalized  so  that  |p4/c|/|E  |  +  ]E  |  =  1  at  x  =  y  =  0. 

x  y 

Broken  Lines  are  for  Negative  Values. 


Figure  19c. 


Constant  Value  Contours  for  Normalized  Field  Component  |Ez/w| 
of  the  TM^  2  Mode  when  h/w  =  3.  The  Fields  are  Normalized  so 
that  jp'V^j/jE  ;  +|E  1^  =  1  at  x  =  y  =  0. 
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Figure  19e.  Constant  Value  Contours  £or  Normalized  Field  Component 

/(Im(p2E JO)2  +  (Im(p2Ev/?))2  -  /(InKp2!!  /se  >)2  +  (Im(p2H  / 
a  y  y  o  a 


x'8e 


of  the  TM.  ~  Mode  when  h/w  -  3.  The  Fields  are  Normalized  so 
■»  ■L»f/ — -n — - r 


that  (p  /c|/lE_r+ |E  |  -  1  at  x-y-0. 

*  y 


x/w 

Figure  19f.  Constant  Value  Contours  £or  Normalized  Field  Component 

|P2/?|/|EJ2+  |Ey|2  •  |p2(se0)_:L|/|Hx|2+|Hy|2  of  the 
Node  when  h/v  ■  3.  The  Fields  are  Normalized  so  that 


SECTION  V 


CONTINUOUS  SPECTRUM  CONTRIBUTION 


In  this  section,  an  asymptotic  estimation  of  the  continuous  spectrum 
contribution  to  the  TM  fields  at  a  given  to  will  be  presented. 

From  equations  2  and  19,  the  continuous  spectrum  contribution  to  E 

z 

for  a  fixed  to  is  given  by 

E^(x,y,z,ju)  |  j*  jjf  J  K0(p/x-x')2+  (y-h)2  )- Ko(p/(x-x')2+  (y+h)2  ) 

n_  -w 

x  f”e(x'/w)e^zdx,dc  (22) 

where  p  *  -  j/?  +  to  /c  ,  B_  is  the  contour  shown  in  figure  3,  and  f  (5) 
is  calculated  from  equations  11  and  18.  Along  B_,  it  can  be  shown  that 
the  p-values  above  and  below  the  branch  cut  are  related  by  p (above)  * 
e^p  (below) .  Thus,  with  the  variable  change  from  t  to  -  ju/c  -  k, 
equation  22  becomes 


E®(x,y,z,ju)  *  e~^“Z^C  j*  dx'  |  die  e  KZf  e(x'/w) 

-w  0 


(23) 


I  (p«4x-x')2+  (y-h)2  )  -  I  (p/(x-  x')2+  (y  +  h)2 


where  IQ  is  the  modified  Bessel  function  of  the  first  kind  and  the  parameter 
p  is  given  by  p  ■  -  j/c2  +  2j<w/c. 

The  asymptotic  behavior  for  large  z  of  the  continuous  spectrum  contri¬ 
bution  Ec(x,y,z,1u)  can  now  be  estimated.  Due  to  the  exponential  term 
exp(-KZ),  it  is  clear  that  the  K-integral  in  equation  23  comes  mainly 
from  the  region  where  1  >.  kz  >.0.  For  field  points  in  the  working  volume 

of  the  simulator,  the  arguments  of  the  Bessel  functions  are  small  for 

2  2  2  2  2 

1  >  «  >  0,  provided  that  z  »  w(h  +w  )/c  and  z  »  (h  +u  ).  The  small- 
argument  expansion  can  be  applied  to  the  Bessel  functions  to  get  (ref.  9) 


(24) 


Io(p/Oc- x’)2  +  (y-h)2  )-  ijp^x-x’)2  +  (y  +  h)2  ) 
«  (k2  +  2Jicw/c)yh 


— ^ 

For  the  term  £  (5).  one  obtains,  under  the  same  conditions,  from  equations 

11  and  18 


T 

where  [T2n(£)]  is  the  transpose  of  the  column  vector  [T^U)]  (i.e. ,  a  row 
vector  with  T2q(0  as  the  elements). 

Combining  equations  23,  24  and  25,  one  has 


E^x.y.z.ju)  a  ^  e~^z,c  l 

n-Q 


i1  +  5no  l£l  S)’1 


[  d5'  — g-  f  die  s“ee"'e*(K2  +  2j icu/c) 

"I  ATT1  0 

2  2  2  2  2 
for  z  >>  ui(h  +w  )/c  and  z  »  h  +  w 


(26) 


Now,  if  one  assumes  that  SQ  has  no  singularities  close  to  the  branch  cut, 
which  is  generally  true,  equation  26  immediately  becomes 


E*(*.y»*,J«> 


in 


4h  \_1  „-e  -juiz/c 

fw  J  ao  e 


(27) 
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From  equation  27,  it  is  seen  that  the  continuous  spectrum  contribution 

_2 

to  E  (x,y,z,ju)  decays  as  z  when  the  wave  propagates  along  the  +z  direc- 
z 

tion.  The  above  asymptotic  estimation  can  also  be  applied  to  other  field 

components  to  obtain  similar  results.  Hence,  within  the  region  where 
2  2222222 

z  »  yhww/c,  (h+w),«a(h+w)/c  ,  the  continuous  spectrum  contribution 
is  negligible  compared  to  the  TEM  mode  contribution. 
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APPENDIX 


MATRIX-EQUATION  FORMULATION  FOR  THE  TE  FIELDS 

In  this  appendix,  the  matrix  equation  for  the  TE  fields  of  the  two 
parallel-plates  will  be  formulated.  The  procedure  will  be  the  same  as 
that  used  for  the  TM  fields  given  in  section  III. 

Under  the  source  condition  6,  the  coupled  differential-integral 
equations  5  are  simplified  to  the  single  equation 


_  p2^j  {  G(x,h,x' ,h;p)  -  G(x,h,x'  ,-h;p)  }  g'^x’Jdx'^ 

-  g'°(x)  for  |  x  |  <.  w  (Al) 

To  solve  equation  Al,  one  first  integrates  the  equation  to  yield  the 
following  integral  equation 

,v 

I  {  G(x,h,x' ,h;p)  -  G(x,h,x' ,-h;p)  }  g"°(x’ )dx' 

-w 

*  B  °sinh(px)  +  (2p)  ^sinh(p|x- x' |  )e  °(x,)dx'  for  |x|  <_w  (A2) 

-w 

where  the  integration  constant  B  0  will  be  determined  from  the  edge  conditions 

Equation  A2  has  the  same  form  as  equation  7.  Hence,  by  following  the 
same  procedure  used  to  solve  equation  7,  one  can  transform  equation  A2  into 
a  Fredholm  integral  equation  of  the  second  kind  given  by 

g“°(5)  +  f  (ku.c’jy)  -  LU.C’JY)) 


for  k  <.  1 


where  the  variables  are  normalized  according  to  equation  8,  K(5,£';y) 
and  L(£,£';y)  are  defined  in  equation  10,  and 

v’°(5)  -  v2  f1  (2Y)"1sinh  (y|5-«'|) 

J-1 

The  unknown  quantities  in  equation  A3  are  g  °(£)  and  B  By 
applying  the  edge  conditions,  one  can  derive  a  relationship  between  them. 
The  edge  and  symmetry  conditions  for  g  °(£)  suggest  the  following  expansion 

g-°«)  -  l  Cm®  (m> 

m“o 

After  using  the  expansion  A4  in  equation  A3  it  is  observed  that  both _ 

left-hand  and  right-hand  sides  of  equation  A3  contain  the  term  £//l-  C2 
which  is  the  dominant  term  when  £  -►  ±1.  By  letting  £  -►  ±1,  one  can  equate 
the  left-hand  and  right-hand  coefficients  of  the  term  £//l-  £*  and  thus 
obtain 


Or,  after  simplification 


9  0  ■  l  a°g  0  -  v  ° 

m:0  A  ed8e 


(A5) 


where 


«  2(irYl0(Y))-1  |  (m'(2h;y)  -  N' (2n;Y))R°(h)dn 

l-|n|  (l- (n' - n)2)u2m+1(n' -  n) 


R°(n) 

a 


l 


-Mn|  A  -  (n*+n)2)(l-  (n'-n)2) 
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*2,,.  ■  f  '''°'(E')(l-5'2)'‘i  dV 

-1 

Having  determined  B  °  in  terms  of  g  the  matrix  equation  for  g  ° 
can  then  be  derived.  By  removing  the  term  C//l-  C  from  both  sides  of 
equation  A3,  it  is  easy  to  see  that  the  equation  can  be  cast  into  the 
following  form 


g-°(o  + 


A-i2  f  {p(c,C;y)  -  Q(5,5,;Y)}g“°(C,)dC’ 

_ -1 

A-  c2  (  b“°h°(5)  +  v"°(o) 


where 


-1  1  “  5 


2l  f1  coshCYg') 

^  1,  / - r2 

_1  («-5')A-5 


dC* 


{p(C,C';y)  -  Q(5.5';y)} 


If1  M'q”-i-’;Y)  -  N(C”-C’; 

-1  (c-5”)^-V'2 


^d£" 


v“°(£) 


-h: 


y"°  (gf) 


'-1 


(€-  5 


di' 


By  inserting  the  expansion  A4  into  equation  A6  and  using  the  orthog' 
onality  of  the  Chebyshev  polynomials,  one  finally  obtains  the  desired 


algebraic  equations  for  gm°  given  by 


g“°  +  l  (  P°  -  Q°  -  H°a°  )g-°  -  -  v“°  H°  +  V"° 
®n  L  \  nm  Tim  n  m  )  m  edge  n  n 


m*o 
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Or,  in  matrix  form 


[d  +  P°  -  Q°  -  H°a°  1  r  g_°  1  -  [-  v“®  H°  +  V"°  1 
L  nm  nm  Tim  n  m  J  L  m  J  (.  edge  n  n  J 


where 


*  C  *  f  £  ("  /iT^U2„+l<»/f75r2  W5'>  « 

{?«,«' rr>  -  Q(£,r  ;y)}drdc 

m  ~  ^  j  (h' (2n;y)  -  N' (2n;y)) F^m(n)dn 


nm 


H°.2[l 

n  ir  J_i 


n-  z2  u2n+1(c)H°(?)d5 


8  f^2 

— ^  cosh(y  cos  4>)cos(  2(n  +  l)<|>)d<j> 

'a 


r°  - 1 

n  it 


rl 


A-z2  u. 


«’  U2n+l(5)V"°(C)d5 


8  f1  -o’ 

"L v 


(C,)T2n+2(5,)(1-5,2)'*S  dr 


On) 

nm 


fl- 1 n |  (l-  (n,-n)2)u2m+1(n,-n)T7n+2(n,+Ti) 


2n+2 


■1+ln*  /( 1  -  (n'  -  n)2 )( l -  (rT+”n)2)~ 


dn' 


n+1  m 

I  l 

k-0  5.-0  nK 


H+k 


j;  ( b  (2k,  2JI+1)  -  nb.(2k,  2£+3) ) 
j-0  x  2  3 


n2t  +  2k-2j  +  lDj(ii)  -  bl+k+1(2k.  2t«)nDt+k+l(n) 


(A8) 
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o  ,  .  v  m  -  5-  (n+Hl)!  .2%  +  1 

um£  =  '  (m-  i) !  (2£  + 1) ! 


o  /  .^n-k+l  (n  +  k)!(n+l) 
nk  ~  1  ’  (n  -  k  +  1) !  (2k) ! 


and  bj(k,A),  (n)  are  given  in  equations  16  and  17;  vedge  and  am  are  given 
in  equation  A5  with 

in 

R°(n)  =  l  U°£  [  (bj(0»  n+1)  -  nVj(0,  2JH-3))n2i-2j  + V(n) 

W°- 


The  solution  of  equation  A8  is  simply  given  by 

'  [  g“°  1  =  \  5  +  P°  -  Q°  -  H°a°  1  \  V“°  -  v~°  H°  1 

J  [  na  om  tu&  n  mJ  L  n  edge  n  J 


This  solution  can  be  used  in  equations  A4,  4  and  2  to  calculate  the  TE .fields 
of  the  two-parallel-plate  simulator. 
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